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We review a recent approach to determine the minimal spatial length scales on which local tem- 
perature exists. After mentioning an experiment where such considerations are of relevance, we first 
discuss the precise definition of the existence of local temperature and its physical relevance. The 
approach to calculate the length scales in question considers homogenous chains of particles with 
nearest neighbor interactions. The entire chain is assumed to be in a thermal equilibrium state and 
it is analyzed when such an equilibrium state at the same time exists for a local part of it. The 
result yields estimates for real materials, the liability of which is discussed in the sequel. We finally 
consider a possibility to detect the existence or non-existence of a local thermal state in experiment. 
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I. INTRODUCTION 

Large systems in an equilibrium state may, despite 
their very large number of degrees of freedom, be char- 
acterized by only very few quantities. For example, an 
ideal gas is described by the simple "thermal equation of 
state" pV = nksT, where p is the pressure of the gas, V 
its volume, n the number of particles it contains, T its 
temperature and ks Boltzmann's constant. In physics 
one refers to this kind of description as a thermodynam- 
ical description. 

How can such an extremely reduced description be jus- 
tified? The reason why a thermodynamical description 
works so well for equilibrium states is that, with increas- 
ing number of particles in a system, a dominant part of 
microstates have the same macroscopic properties. Mi- 
crostate refer here to a description where all degrees of 
freedom are specified. As a result, thermodynamical be- 
havior becomes "typical" . 

In a more mathematical language this fact is called 
the existence of the Thermodynamic Limit, which merely 
means that intensive quantities such as the energy per 
particle approach a limiting value that does no longer 
depend on the detailed configuration of the system as 
its size increases. For example the energy per particle 
of a very large piece of solid does no longer depend on 
whether this piece is lying on a table or is immersed in 
a bucket full of water, provided it is in an equilibrium 
state, i.e. has the same temperature as its surrounding. 

Obviously, in the Thermodynamic Limit, the differ- 
ence between particles inside the solid and those on the 
surface, which interact with the surrounding, becomes 
negligible. That is why the size of the solid is important. 
Let us assume the solid had the shape of a sphere and the 
density of the particles was uniform within it. Since the 
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surface of a sphere with radius r is Anr 2 and its volume 
is (47r/3)r 3 , the ratio of particles sitting on the surface 
over the total number scales as 1/r and thus becomes 
negligible as r goes to infinity. Such a type of scaling 
does apply not only to spheres but also to more general 
geometries. 

To analyze the existence of local temperatures, small 
parts of large systems are of interest. These parts do in- 
evitably interact with their surrounding. For short range 
interactions between the constituent particles, again, 
only those particles sitting on the boundary of the con- 
sidered part interact with the environment. Thus, the 
described scaling properties immediately give rise to the 
following question: How large do parts of those systems 
have to be to permit a local thermodynamical descrip- 
tion, i.e. a thermodynamical description of the part 
alone? 

For a long time, the problem, besides being funda- 
mental, may have been of purely academic interest, 
since thermodynamics was only used to describe macro- 
scopic systems, where deviations form the Thermody- 
namic Limit may safely be neglected. However, with the 
advent of nanotcchnology, the microscopic limit of the 
applicability of thermodynamics became relevant for the 
interpretation of experiments and may in the near future 
even have technological importance. 

In recent years, amazing progress in the synthesis and 
processing of materials with structures on nanometer 
length scales has been made 0, 0, IE 0- Experimen- 
tal techniques have improved to such an extent that the 
measurement of thermodynamic quantities like temper- 
ature with a spatial resolution on the nanometer scale 
seems within reach [1110. To provide a basis for the 
interpretation of present day and future experiments in 
nanoscale physics and technology and to obtain a bet- 
ter understanding of the limits of thermodynamics, it is 
thus indispensable to clarify the applicability of thermo- 
dynamical concepts on small length scales starting from 
the most fundamental theory at hand, i. e. quantum 



mechanics. In this context, one question appears to be 
particularly important and interesting: Can temperature 
be meaningfully defined on nanometer length scales? 

Why should we care about the non-existence of lo- 
cal temperature? There are at least three situations for 
which this possibility needs special attention: One obvi- 
ous scenario refers to the limit of spatial resolution on 
which a temperature profile could be defined. However a 
spatially varying temperature calls for non-equilibrium - 
a complication which we will exclude here. A second ap- 
plication concerns partitions on the nanoscale: If a mod- 
ular system in thermal equilibrium is partitioned into two 
pieces, say, the two pieces need no longer be in a canonical 
state, let alone have the same local temperature. Finally, 
local physical properties may show different behavior de- 
pending on whether the local state is thermal or not. 

The existence of thermodynamical quantities, i.e. the 
existence of the Thermodynamic Limit strongly depends 
on the correlations between the considered parts of a sys- 
tem. As mentioned above, with increasing diameter, the 
volume of a region in space grows faster than its surface. 
Thus effective interactions between two regions, provided 
they are short ranged, become less relevant as the sizes 
of the regions increase. This scaling behavior is used 
to show that correlations between a region and its envi- 
ronment become negligible in the limit of infinite region 
size and that therefore the Thermodynamic Limit exists 
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To explore the minimal region size needed for the appli- 
cation of thermodynamical concepts, situations far away 
from the Thermodynamic Limit should be analyzed. On 
the other hand, effective correlations between the consid- 
ered parts need to be small enough 0, 0] . 

The scaling of interactions between parts of a system 
compared to the energy contained in the parts themselves 
thus sets a minimal length scale on which correlations 
are still small enough to permit the definition of local 
temperatures. Here we review an approach to study this 
connection quantitatively [l3l. Ilif . 



II. MOTIVATION: A THERMAL NANOSCALE 
EXPERIMENT 

In recent years, there has been substantial progress in 
the fabrication and operating of material with structure 
on nanoscopic scales and nanoscale devices. In this con- 
text, several experiments, that study thermal properties, 
have been done. We describe here, as an example, one 
experiment that nicely shows where the existence or non- 
existence of local temperature becomes relevant 0] . 

The experiment studies heat conduction across a car- 
bon nanotube. A sketch of the setup is given in figure ^ 
Two, otherwise thermally well isolated islands are con- 
nected through a carbon nanotube of a few length. 
One island is heated by an electric current that runs 
through a coil with the resistance This island is 

thus at a "hot" temperature Xh. Heat can flow across 
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FIG. 1: Setup of the experiment. Two, otherwise thermally 
well isolated islands are connected by a carbon nanotube. The 
left island is heated by an electric current running through the 
coil with resistance Rh and thus maintained at the tempera- 
ture Th- The temperature of the right island is measured via 
the temperature dependent resistance R s . 




FIG. 2: Picture of the setup. The heated island is in the lower 
left corner and the island where the temperature is measured 
in the higher right corner. Both are connected by a single 
carbon nanotube (With permission from A. Majumdar, Na- 
noengineering Laboratory, University of California, Berkeley) . 

the nanotube to the other island, which is at a lower 
temperature T s . This temperature in turn is measured 
by another coil, the resistance of which R s , depends on 
temperature. Figure [21 shows a picture of this setup. 

At what point is the existence or non-existence of local 
temperatures of relevance for the interpretation of this 
experiment? We know that there is an electric current in 
the coil of the heated island. This current constantly de- 
livers thermal energy to the island. This energy is trans- 
ported across the nanotube to the other island, where we 
observe, that the temperature T s rises. We thus know, 
that the nanotube connects a hot spot Th to a cold spot 
T s . This directly gives rise to the following questions: 
How hot is the nanotube in between? Can we mean- 
ingfully talk at all about temperature of any part of the 
nanotube? The answer to these two questions would clar- 
ify whether and in what sense a temperature profile (sec 
figure |3) could exist for the present setup. 

While a temperature profile can obviously be defined 
and measured in a macroscopic version of the present 
experiment, say two buckets of water at different tem- 
peratures and connected via an iron bar, its existence, 
possible resolution and measurability are completely un- 
clear for the nanoscopic version. 
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FIG. 3: The question, whether a temperature profile exists 
for the nanotube in the present setup, is not clarified. 



To address this question, we first discuss how to define 
the existence of local temperature in the next section. 



III. WHAT IS TEMPERATURE? 

Temperature is one of the central quantities in thermo- 
dynamics and Statistical Mechanics. Let us note here, 
that it is not a direct observable; it is not represented 
by an operator in quantum mechanics (cf. sect. IVIII(I . 
There exist two standard ways to define it: 



A. Definition in Thermodynamics 

The thermodynamical definition is purely empirical. 
Thermodynamics itself is an empirical theory on systems 
whose macroscopic physics can be sufficiently character- 
ized by a set of a few variables like volume, energy and 
the number of particles for example. The values of this 
variables are called a macro state [l^. A system is said 
to be in equilibrium if its macro state is stationary for 
given constraints. As a consequence of this definition, 
a equilibrium state depends on the applied constraints, 
e.g. whether the volume or the energy is kept constant 
etc. An important, special case of constraints, is a bipar- 
tite (or multipartite) system with a fixed total energy, 
where the parts may exchange energy among themselves. 
The parts are then said to be in thermal equilibrium. In 
thermodynamics, temperature is defined by the following 
property: 

Definition: Two systems that can exchange energy 
and are in thermal equilibrium, have the 
same temperature. 

To fix a temperature scale, a reference system is 
needed. The simplest choice for this reference system 
is the ideal gas (cf. sec. HJ, where temperature may be 
defined by 



T = 



P V 
n fee 



Of course the above definition is unambiguous only if 
the states of thermal equilibrium form a one-dimensional 
manifold ^||. Only then, a single parameter is sufficient 
for their characterization. This parameter is the temper- 
ature T. 



B. Definition in Statistical Mechanics 

In statistical mechanics, temperature is defined via the 
derivative of the entropy S with respect to the internal 
energy E. In quantum mechanics the entropy can be 
defined according to von Neumann as 



S = — /cb Trplnp , 



(2) 



it is a measure of the amount of possible pure states, the 
system could be in. With this definition, entropy always 
exists, but shows its standard properties, e.g. extensivity, 
only in the Thermodynamic Limit 0. 

In statistical mechanics, an equilibrium state is defined 
to be the state with the maximal entropy, that is the state 
with the maximal amount of accessible pure states. 

For systems that interact with their surrounding, such 
that they can exchange energy with it but have a fixed 
expectation value for the energy, the equilibrium state is 
a so called canonical state, described by a density matrix 
of the form 



P 



exp{-f3H) 



(3) 



where the partition sum Z normalizes p such that Tr p = 
1. 

In quantum mechanics the internal energy is given by 
the expectation value of the energy, 



E = Tr pH , 



(4) 



where the Hamiltonian is the energy operator of the iso- 
lated system at hand. It does not contain any interac- 
tions of the system with its environment. The internal 
energy is therefore a property of the system itself, it only 
depends on the state of the system and not on the state 
of the environment. 

Temperature is then defined by 



1 

T 



dS_ 



(5) 



(1) 



which in turn exists as long as the entropy S is a func- 
tion of the internal energy E. However, the notion of 
temperature, as defined in equation (J5J, just like entropy 
shows its characteristic thermodynamical properties (see 
above) only for equilibrium states |sL Il6l lri| . 



C. Local Temperature 

Local temperature is, by definition, the temperature of 
a part of a larger system. Hence, this subsystem is not 
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isolated but can exchange energy with its surrounding. 
On the other hand we limit our considerations to cases 
without particle exchange. Hence, the following conven- 
tion appears to be reasonable: 

Definition: Local temperature exists if the consid- 
ered part of the system is in a canonical 
state. 

Note: While a local state can always uniquely be de- 
fined by tracing out the rest of the system, this definition 
calls, in addition, for the (approximate) existence of some 
local spectrum. 

Besides being based on statistical mechanics there are 
further practical reasons for this definition: The canoni- 
cal distribution is an exponentially decaying function of 
energy characterized by one single parameter, tempera- 
ture. This implies that there is a one to one mapping 
between temperature and the expectation values of ob- 
servables, by which it is usually measured. Tempera- 
ture measurements via different observables thus yield 
the same result, contrary to distributions with several 
parameters. 

This is a basic property of systems that can be charac- 
terized by thermodynamic description. The temperature, 
if it exists, describes a system in a sufficiently complete 
way, such that several properties of it can be predicted if 
one only knows its temperature (see also sec. IVffffl . 

Why does the distribution need to be exponentially 
decaying? In large systems with a modular structure, 
the density of states is a strongly growing function of en- 
ergy 01 • The product of the density of states times an 
exponentially decaying distribution of occupation proba- 
bilities will thus form a strongly pronounced peak at the 
internal energy E, see figure 0] 
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FIG. 4: The product of the density of states n(E) times the 
occupation probabilities (</?| p \tp) forms a strongly pronounced 
peak at E = E. 

If the distribution was not exponentially decaying, the 
product of the density of states times the distribution 
would not have a pronounced peak and thus physical 
quantities like energy could not have "sharp" values. 



IV. GENERAL THEORY FOR THE EXISTENCE 
OF LOCAL TEMPERATURE 



After having introduced and discussed the conception, 
when temperature is defined to exist locally, we now turn 
to describe the approach to analyze its existence. 

Since temperature is defined to exist locally, i.e. for 
a given part of the system we consider, if the respective 
part is in a thermal equilibrium state, it is defined to exist 
on a certain length scale, if all possible partitions of the 
corresponding size are simultaneously in an equilibrium 
state. 

This requirement for the local equilibrium states in 
the parts to exist at the same time needs some further 
discussion: For a given temperature profile, it should not 
make a difference whether the profile is scanned by one 
single thermometer, which is moved in small steps across 
the sample, or whether the profile is measured by several 
thermometers simultaneously, which are located at small 
distances to each other. 

For systems which are globally in a non-equilibrium 
state it is very difficult to decide under what conditions 
equilibrium states show up locally ^g] and only very few 
exact results are known Nonetheless, whenever local 
equilibrium exists, the macroscopic temperature gradient 
is small (ST/T <C I). Here, we restrict ourselves to sys- 
tems which are in a global equilibrium state In these 
situations, subunits of the total system are in an equi- 
librium state whenever their effective interaction is weak 
enough and correlations between them are small so that 
the global thermal state approximately factorizes into a 
product of local thermal states. 

Whenever the macroscopic temperature gradient is 
small (5T/T <C I), one would expect the results to be 
applicable even for situations with only local equilibrium 
but non-equilibrium on the global scale. 

To explore how local temperature can exist, that is 
how small the respective part may be, one needs look at 
parts of different sizes. The idea behind this approach is 
the scaling behavior which ensures the existence of the 
Thermodynamic Limit (cf. section HJ 0,^3 : 

We consider systems that are composed of elementary 
subsystems with short range interaction, for simplicity 
say nearest neighbor interaction. If then n adjoining sub- 
systems form a part, the energy of the part is n times 
the average energy per subsystem and is thus expected 
to grow as the size of the part, n. Since the subsystems 
only interact with their nearest neighbors, two adjacent 
parts interact via the two subsystems at the respective 
boundaries, only. As a consequence, the effective cou- 
pling between two parts is independent of the part size n 
and thus becomes less relevant compared to the energy 
contained in the parts as their size increases. 
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FIG. 5: Groups of n adjoining subsystems are formed. 



A. The Model 

As models we consider here homogeneous (i.e. trans- 
lation invariant) systems with nearest neighbor interac- 
tions which we divide into identical parts. The Hamilto- 
nian of the system thus reads 



(6) 



where the index i labels the elementary subsystems. Hi 
is the Hamiltonian of subsystem i, ij^+i the interaction 
between subsystem i and i + 1 and periodic boundary 
conditions are assumed. 

Now Nq groups of n subsystems each (index i — > (/i — 
l)n + j ; jji = 1, . . . , Nq ; j = 1, . . . , n) are formed and the 
Hamiltonian is split into two parts, 



H = H +I, 



(7) 



where Ho is the sum of the Hamiltonians of the isolated 
groups, 



with 



(8) 



n—1 



and / contains the interaction terms of each group with 
its neighbor group, 



N a 

I — ^ ^ I/j.n,tin-\-l 



(9) 



The eigenstates \a) of the Hamiltonian Ho, 
Hq | a) = E a \a), are products of group eigenstates 
of the individual groups, 



Nq 



with W M |a M ) = £ M |a M ) , (10) 



where is the energy of one subgroup only and 

B. Thermal State in the Product Basis 

To test whether a part H^ is in a thermal state, we 
have to calculate its reduced density matrix by tracing 



out the rest of the system. This trace can only be per- 
formed in the basis formed by the states \a), ifTUI) . We 
thus have to write the global equilibrium state © in this 
basis. Denoting the eigenstates and eigenenergies of the 
global Hamiltonian with Greek indices, \ip), \ip) and E, 
E^p, the global equilibrium state p reads 



(ii) 



in the global eigenbasis and the diagonal elements in the 
product basis are 



(a\p\a) 



e -f3E 

Wa (E)—dE, (12) 



where the original sum has been replaced by an in- 
tegral over the energy. 

w a (E) is the probability to obtain an energy value be- 
tween E and E + AE if the total energy H is measured 
for as system in the state \a), i.e. 



(13) 



{\ip):E<E v <E+AE} 



where the sum runs over all states \tp) with energy eigen- 
values E v in the respective energy range and AE is small. 
Eq is the energy of the ground state and E\ the upper 
limit of the spectrum, which should be taken to be infi- 
nite if the spectrum does not have an upper bound. 

We thus have to know the distributions w a {E) in order 
to be able to compute the reduced density matrices of the 
groups and to test whether they are of canonical form. 

Fortunately, one can indeed show that there exists a 
quantum central limit theorem for many particle systems 
with nearest neighbor interactions |20l |21| . Therefore, 
in the limit of infinitely many groups, w a takes on the 
following form: 



lim 



w a (E) 



E 



/ 2tt A a 



exp 



E a r 



(14) 



where E a , the expectation value of H in the state |o), 
and A^, its variance, read 



E a = (a\H\a) 



A 



(a\H z \a) 



and 
(a\H\a) 



(15) 
(16) 



Note here, that the limit of infinite number of groups 
is taken while the size of each individual group remains 
finite. 

For the theorem to hold, two further conditions have 
to be met: The energy of each group including its inter- 
actions with the neighboring group has to be bounded 
and the variance A a has to grow faster than NqC for 
some positive constant C . In scenarios, where the en- 
ergy spectrum of each elementary subsystem has an up- 
per limit, such as spins, the first condition is met a priori. 
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For subsystems with an infinite energy spectrum, such as 
harmonic oscillators, the present analysis is restricted to 
states where the energy of every group, including the in- 
teractions with its neighbor groups, is bounded. Thus, 
the considerations do not apply to product states \a), 
for which all the energy was located in only one group 
or only a small number of groups. The number of such 
states is vanishingly small compared to the number of all 
product states. 

The expectation value of the entire Hamiltonian H in 
the state |a), E a , is the sum of the energy eigenvalue 
of the isolated groups E a and a term that contains the 
interactions, 



(17) 



Therefore, the two quantities e a and A 2 can also be ex- 
pressed in terms of the interaction (see eq. J7J)) only, 



= (all la) 



and 



A„ = (a\I 2 \a)-(a\I\a) 



(18) 
(19) 



meaning that e a is the expectation value and A 2 the 
squared width of the interactions in the state \a). Note 
that e a has a classical counterpart while A 2 is purely 
quantum mechanical. It appears because the commuta- 
tor [H , Hq] is nonzero, and the distribution w a (E) there- 
fore has nonzero width. 

Applying equation l|14|) to calculate the integral in 
equation {HJ yields for N G > 1, 

1 / /3 2 A 2N 
(o| p\a) = — exp -p (E a + e B ) + ^— - I 



Z 
1 



t M,l E °- Ea Z £a + PAr )- (20) 



erfc 



V2A a 

E x - E a -e a + (3A 2 a 
V2A a 



where eifc(x) is the conjugate Gaussian error function 
, The second error function appears only if the energy 
is bounded and the integration extends from the energy 
of the ground state Eq to the upper limit of the spectrum 
Ex. 

Note that the arguments of the conjugate error func- 
tions grow proportional to s/Nq or stronger, therefore 
the asymptotic expansion of the latter |22| may be used 
for N G > 1. 

The off diagonal elements (a\p\b) vanish for 
\E a — E \ > A a +Af, because the overlap of the two Gaus- 
sian distributions becomes negligible. For \E a — Eb\ < 
A a + Ab, the transformation involves an integral over 
frequencies and thus these terms are significantly smaller 
than the entries on the diagonal. 



C. Conditions for Local Thermal States 

We now test under what conditions the density ma- 
trix p may be approximated by a product of canoni- 
cal density matrices with temperature (3i oc for each sub- 
group p = 1,2,...,Nq. Since the trace of a matrix 
is invariant under basis transformations, it is sufficient 
to verify the correct energy dependence of the product 
density matrix. If we assume periodic boundary con- 
ditions, all reduced density matrices are equal and if 
they were canonical their product would be of the form 
(a\ p \a) oc exp(— j3\ oc E a ). We thus have to verify whether 
the logarithm of rhs of equation H20|) is a linear function 
of the energy E a defined in equation (|l(Jfl and below, 



In ((a| p \a)) w -0\ oc E a + c, 



(21) 



where (3\ oc and c are constants. 

Applying the asymptotic expansion of the conjugate 
error function to l|20[l shows that equation l|21l) can only 
be true for 



E a 



En 



> 



— A 2 

2 a 




and 



(22) 



(23) 



where c\ and ci are constants. These two conditions 
constitute the general result of this section. 

Note that e a and A 2 need not be functions of E a and 
therefore in general cannot be expanded in a Taylor se- 
ries. 

Temperature becomes intensive, if the constant C\ van- 
ishes, 



Icil « 1 



Aoc = ■ 



(24) 



Even if this was not the case, temperature might still 
exist locally. 

For the existence of local temperature, one should only 
require that the diagonal elements (|12f> are canonically 
distributed in an appropriate energy range, E m { n < E a < 
£ max . As described in the previous section, the den- 
sity of states r)(E) is, for large modular systems, an ex- 
ponentially growing function of energy and its product 
with the exponentially decaying canonical distribution 
{ip \ p \ (p) forms a strongly pronounced peak at the expec- 
tation value of the global energy E, see figure 0J 

If the diagonal elements l|12|l are canonically dis- 
tributed in an energy range, that is centered at this peak 
and is large enough to entirely cover it, all observables 
with non-vanishing matrix elements in that range show 
the same behavior as for a true canonical distribution. 
Observables which are not of that kind are in general not 
of interest. If one considers for example 1 kg of iron at 
300 Kelvin with an average energy of roughly 130 kJ, one 
is usually not interested in processes, that take place at 
energies of 0.1 kJ or 10 5 kJ. 
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Therefore a pertinent and "safe" choice for the energy 
range E min < E a < E max should be 

E min = max ( [E a ] ■ , -E + E ) 

V a J (25) 

^max = niin ([£ a ] max , otE + E ) , 

where a>l and E will in general depend on the global 
temperature (3. In equation [-E>] min and [-E>] max 

denote the minimal and maximal values E^ can take on. 

In ({a\ p \a)) 




E 



FIG. 6: \n((a\p\a)) for p as in equation 1201 (solid line) and a 
canonical density matrix p (dashed line) for a harmonic chain. 

Figure |H1 shows the logarithm of equation (|20|l and the 
logarithm of a canonical distribution with the same (3 for 
the example of a harmonic chain. The actual density ma- 
trix is more mixed than the canonical one. In the interval 
between the two vertical lines, both criteria (1221 and 12: il) 
are satisfied. For E < E\ ow <|22|) is violated and (|23fl for 
E > -Ehigh- To allow for a description by means of canon- 
ical density matrices, the group size needs to be chosen 
such that E iow < E min and E high > -E m ax- 

For a model of the class considered here, the two condi- 
tions l|22ll and (|23[) must both be satisfied. In the follow- 
ing sections, these fundamental criteria will be applied to 
a concrete model. 



V. HARMONIC CHAIN 

We consider a harmonic chain of Nn ■ n particles of 
mass m and spring constant ^Jtuloq [jjj,|23]. In this case, 
the respective terms in the Hamiltonian JfJJl read 

Hi = y^ + y^k 2 (26) 
h,i+l = -mujgqiq i+ i , (27) 

where pi is the momentum of the particle at site i and 
qi the displacement from its equilibrium position i ■ ciq 
with ao being the distance between neighboring parti- 
cles at equilibrium. We divide the chain into Nq groups 
of n particles each and thus get a partition of the type 
considered in section llVl 



The Hamiltonian of each group is diagonalized by a 
Fourier transform and the definition of creation and an- 
nihilation operators a\. and ak for the Fourier modes • 
In this way we get 

n g / 1\ 

E a = with Ep = £>fc I ngfa) + - J , (28) 

fj.=l k ^ ' 

where k = nl/(ao (n +1)) (/ = 1,2,..., n) and the fre- 
quencies cok are given by io\ = 4w 2 sin 2 (kao/2). n^.(fi) 
is the occupation number of mode k of group (i in the 
state | a). We choose units, where H = 1. Let us first 
verify, that the central limit theorem (|14|l applies to this 
model, i.e. that the required assumptions (given below 
equation l(Trj|) ') are met. 

To see that A 2 grows faster than NqC (C > 0), one 
needs to express the group interaction V {q^m q_^n+i) in 
terms of a\ and a^; this yields A 2 = Y^lii^tfi where 
A 2 > 0, implying that the assumption is met. 

Since the spectrum of every single oscillator is infi- 
nite, the requirement that the energy per group should 
be bounded can only be satisfied for states, for which 
the energy of the system is distributed among a relevant 
fraction of the groups. As discussed in section HV1 states 
where this is not the case constitute only a negligible 
fraction of all product states \a). 

The expectation values of the group interactions 
fea. (|18|l % ) vanish, e a — 0, while the widths A 2 accord- 
ing to equation depend on the occupation numbers 
rik(n) and therefore on the energies E^. To analyze con- 
ditions (|22|l and (|23|l . one makes use of the continuum 
or Debye approximation |24j . requiring n 1, ao <C /, 
where I = noo, and the length of the chain to be finite. 
As will become clear below, the resulting minimal group 
sizes n m in are larger than 10 3 for all temperatures and the 
application of the Debye approximation is well justified. 

Using this approximation we now have Uk — vk with 
the constant velocity of sound v = loq ao and the width 
of the group interaction reads 

A 2 I = 4^^+i- (29) 

where n + 1 sa n has been used. 

The relevant energy scale is introduced by the thermal 
expectation value of the entire chain 

_ /T x2 f e/T 

E = E + N G nk B ei^-j J ^—\ dx ' ( 3 °) 
and the ground state energy Eq is given by 

£o = ^e(|) 2 jf e/T |^ = «, (si) 

where is the Debye temperature [24|. 

Inserting equation l|3U|l and (|31|) into equations 1221) 
and (|23|l , taking into account l|25|) , one can now calculate 
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the minimal n for given a, and T. In doing so, one 
needs to introduce another accuracy parameter 5 which, 
for the rhs of eq. (|23|) , quantifies how much smaller terms 
quadratic and higher order in E a are compared to the 
zero order and linear ones. More precisely, 6 is the ratio 
of the higher order terms to the (at most) linear ones. 

Figure [7| shows n m i n for a = 10 and S = 0.01 as a 
function of T/Q. 
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FIG. 7: in in ELS £1 function of T/Q for a harmonic chain as 
determined by eq. 1221 (solid line) and as determined by eq. 
\2'M (dashed line) for a = 10 and 5 = 0.01. Local temperature 
exists in the shaded region. 

For high (low) temperatures n mul can be estimated by 



^ for T > 9 



(32) 



In addition, local temperatures are equal to the global 
one whenever they exist, f3\ oc = (3, implying that tem- 
perature is intensive (see eq. JUJl). 

In the following section the results obtained above will 
be applied to real materials. 



VI. ESTIMATES FOR REAL MATERIALS 

Thermal properties of insulating solids can successfully 
be described by harmonic lattice models. Probably the 
best known example of such a successful modeling is the 
correct prediction of the temperature dependence of the 
specific heat based on the Debye theory |2^|. Therefore 
one would expect the present approach to give reasonable 
estimates for real materials, too. 

We thus take the results obtained in section for the 
harmonic chain and insert the corresponding parameters, 
in particular the Debye temperature which can be found 
tabulated |24j . One obtains a length scale by multiplying 
n m ; n with the corresponding lattice constant. The min- 
imal length scale on which intensive temperatures exist 
in insulating solids should thus given by 



where ao is the lattice constant, the distance between 
neighboring atoms. 

Since n m ; n has been calculated for a one dimensional 
model the results we obtain here should be valid for one 
dimensional or at least quasi one dimensional structures 
of the respective materials. Let us consider two examples: 

Silicon is used in many branches of technology. In 
its crystalline form, it has a Debye temperature of f=a 
645 K and its lattice constant is ao ~ 2.4 A. Using these 
parameters, figure |S] shows the minimal length-scale on 
which temperature can exist in a one-dimensional silicon 
wire as a function of global temperature. Here, the accu- 
racy parameters a (see eq.(J2SJ) and S (see below ea.(|31|)l 
are chosen to be a — 10 and 8 = 0.01. Local temperature 
exists in the shaded area. 
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FIG. 8: Z m in as a function of temperature T for crystalline 
silicon, ao ~ 2.4 A, ~ 645 K, accuracy parameters a = 10 
and 5 — 0.01. Local temperature exists in the shaded area. 



Recently, carbon has been investigated for the fabrica- 
tion of nano-structured devices [25ll26|. In particular, we 
consider carbon nanotubes here, which are widely used in 
nano-technological experiments. Carbon nanotubes have 
diameters of only a few nanometers. Measurements of 
their specific heat have shown, that their thermal prop- 
erties can be accurately modeled with one-dimensional 
harmonic chains |27| . The presented results can thus be 
expected to be accurately applicable to them. Carbon 
nanotubes have a Debye temperature of w 1100 K and 
a lattice constant of ao ~ 1.4 A. 

Figure|5|shows the minimal length-scale on which tem- 
perature can exist in a carbon nanotube as a function 
of global temperature. It provides a good estimate of 
the maximal accuracy, with which temperature profiles 
in such tubes can be meaningfully discussed 0. Again, 
the accuracy parameters a and 6 are chosen to be a = 10 
and S = 0.01. Local temperature exists in the shaded 
area. 

Of course the validity of the harmonic lattice model 
will eventually break down at high but finite temper- 
atures. The estimates drawn from the considered ap- 
proach, in particular the results presented in figures 
and |3 will then no longer apply. 
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FIG. 9: / m in as a function of temperature T for a carbon 
nanotube. ao ~ 1.4 A, ~ HOOK, accuracy parameters 
a = 10 and 5 = 0.01. Local temperature exists in the shaded 
area. 



VII. DISCUSSION OF THE LENGTH SCALE 
RESULTS 

The length scales one obtains here are, in particular 
for low temperatures, surprisingly large. One might thus 
wonder whether the approach really captures the rele- 
vant physics. Let us therefore discuss some possible lim- 
itations: 

Firstly, one may argue that taking the limit of an in- 
finite number of groups, as required for the central limit 
theorem, will not correspond to the physically relevant 
situations. However, having in mind that we intended to 
analyze when a small part of a larger system can be in a 
thermal state, taking this limit should be well justified. 

Secondly, only one-dimensional models were consid- 
ered. A real physical system, even if it is of a very prolate 
shape, is always three-dimensional. A generalization of 
the approach to those models is thus of high interest. 
However let us stress here, that the general conditions 
(|22|l and ll'.'il) apply to systems of arbitrary dimension, it 
is only the application to specific models, which needs to 
be generalized. 

Furthermore, the harmonic chain is an exactly diag- 
onalizable model, this means that no phonon scattering 
does occur. The purely harmonic model does, for exam- 
ple, not predict any expansion or shrinking of the mate- 
rial caused by heating or cooling. It is therefore possible 
that the harmonic model may fail to give reliable results 
for our present investigation. In particular at low tem- 
peratures, entanglement [2^, |2^, plays an important 
role and this effect can be highly non-linear. 

Finally, one might speculate whether the length scales 
could significantly change if the assumption of a global 
equilibrium state was relaxed. This possibility of course 
exists, nonetheless one would expect the estimates to still 
apply as long as temperature gradients are small. Imag- 
ine, there are two baths attached to the ends of the con- 
sidered harmonic chain of section [VJ If both baths have 
the same temperature, the chain is in a "global" equilib- 
rium state and the present results are applicable. If one 
now continously increased the temperature of one bath, 



the density matrix of the chain would change continously, 
too. Hence, the minimal group sizes would also change 
continuously and the present results should still be good 
estimates, at least for small temperature gradients. 

To clarify whether the above findings are in agreement 
or in conflict with experiments, their measurability needs 
to be considered in more detail. We proceed to do this 
in the next section. 



VIII. CONSEQUENCES FOR MEASUREMENTS 

In this section, we give some examples of possible ex- 
perimental consequences of the local breakdown of the 
temperature concept at small length scales, i.e. of the 
fact that the respective individual subsystems or even 
subgroups of those do not reach a canonical state. 

A. Standard Temperature Measurements 

Temperature is always measured indirectly via observ- 
ables, which, in quantum mechanics, are represented by 
hermitian operators. Usually, one is interested in mea- 
suring the temperature of a system in a stationary state. 
The chosen observable should therefore be a conserved 
quantity, i.e. its operator should commute with the 
Hamiltonian of the system. 

A conventional technique, e.g., is to bring the piece of 
matter, the temperature T of which is to be measured, 
in thermal contact with a box of an ideal gas and to 
measure the pressure p of the gas, which is related to its 
temperature by uk^T = pV (cf. sect. [|J. Since the 
gas is in thermal equilibrium with the considered piece 
of matter, both substances have the same temperature. 
A measurement of p for constant V allows to infer the 
global temperature T of the piece of matter. 

One might wonder, whether a small (possibly even 
nanoscopic) thermometer 5] , which is locally coupled to 
one subsystem of the large chain considered in section llVl 
is capable of measuring a local temperature or whether 
the measurement would show any indications of a possi- 
ble local breakdown of temperature. 

A prerequisite for the above gas thermometer to work 
properly is that the thermometer does not significantly 
perturb the system. For our class of models this means 
that the thermometer system should only be weakly cou- 
pled to the respective subsystem of the chain and that 
it should be significantly smaller than the latter. These 
two requirements ensure that the energy exchange be- 
tween system and thermometer would not significantly 
alter the energy contained in the system. Therefore, this 
measurement scenario can be accurately modeled as fol- 
lows: 

Let the thermometer be represented by a single spin, 
which is locally coupled to a harmonic chain, say. Since 
the coupling is assumed to be weak and the chain is as- 
sumed to be very large and in a thermal state, the present 
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scenario can accurately be modeled with a master equa- 
tion approach [3l|. However, it is a well known result 
of such system bath models, that the reduced density 
matrix of the spin relaxes into a canonical state with the 
temperature being equal to the global temperature of the 
harmonic chain (i.e. the bath). The spin (thermometer) 
thus measures the global temperature of the total chain, 
even for perfectly local coupling. 

As long as the chain is in a global equilibrium state, a 
temperature measurement of this type thus does not have 
any spatial resolution at all. It is only capable of mea- 
suring the global temperature of the chain. Neither can 
any local temperatures be measured nor any signatures 
of their breakdown be detected. 

This conclusion obviously cannot hold anymore for sce- 
narios with only local but no global equilibrium. Macro- 
scopic temperature profiles are routinely measured, with 
the standard technique described above. Whether such 
measurements of temperature profiles are still possible 
for much smaller systems and what their maximally pos- 
sible spatial resolution is in that case should be subject 
of further investigations. 

According to the above considerations, one might think 
that the question of local temperatures for systems in 
global equilibrium was an irrelevant issue since it has no 
observable consequences. This, however is not the case. 
In the following we turn to discuss an example of such 
measurable consequences of the local breakdown of the 
concept of temperature. 



B. Non-thermal Local Properties 

We now turn to observables of the object (chain) itself, 
which could be used to infer local temperatures T\ oc , i.e. 
temperatures of subsystems, provided the subsystems are 
in a canonical state. On the other hand, if the respective 
subsystems are not in a canonical state, this fact should 
then modify the measurement results for those observ- 
ables. 

The minimal group sizes calculated in section II VI de- 
pend on the global temperature and on the strength of 
the interactions between neighboring subsystems. Fur- 
thermore, local temperatures can even exist for single 
subsystems if these are finite dimensional. In the limit- 
ing case of infinite temperature, the density matrix of a 
chain of finite dimensional subsystems is proportional to 
the identity matrix and thus has the same form in every 
basis including the product basis, which in turn implies 
that local temperature would then exist for single sub- 
systems |l lj . For systems composed of finite dimensional 
subsystems, local temperatures do thus exist for single 
subsystems at relatively low global temperatures if the 
coupling is weak, while they do not if the coupling is 
strong. 

Pertinent systems, for which such effects could easily 
by studied, are magnetic materials |3^ . These can in 
many cases be described by spin lattice or spin chain 



models. Since, as we will see below, properties of single 
spins can be infered from measurements of even macro- 
scopic magnetic observables, those materials thus allow 
to study the existence of temperature, as defined by the 
existence of a canonical state, on the most local scale 
possible, i.e. for single spins. 

For a spin- 1/2 system, it is always possible to assign 
a Boltzmann factor and thus a local temperature to the 
ratio of the occupation probability of the higher and lower 
level. Here we consider a homogeneous chain of spin-1 
particles interacting with their nearest neighbors. For 
the interactions, one assumes a Heisenberg model. The 
Hamiltonian of this system reads |33j |: 

n n 

" ••'!>;<•■ • g >1- ( 34 ) 

i=i 3=1 

where <rj, a v - and <r| are the spin-1 matrices. B is an 
applied magnetic field, J the coupling and n the number 
of spins. The coupling J is taken to be positive, J > 
0. The spins thus tend to align anti-parallelly and the 
material is anti-ferromagnetic. The local Hamiltonian 
of subsystem j is Hj = B erj . The system has periodic 
boundary conditions and is thus translation invariant. As 
in the previous sections, the entire system (|34l) is assumed 
to be in a thermal state (see equation ©). 

As an example of an experiment, we will now consider 
two different magnetic observables of a spin-1 system 
with the Hamiltonian l|34fl . The first observable is the 
magnetization in the direction of the applied field, m z , 
which is here defined to be the total magnetic moment 
per particle: 

™* = -(£^)> ( 35 ) 
i=i 

where (O) is the expectation value of the operator O, i.e. 
(O) = Tr(pO). In the translation invariant state p, the 
reduced density matrices of all individual spins are equal, 
and the magnetization (|35fl can be written as 

m z = (a%) , (36) 

for any k = 1,2, ...,n. The magnetization, although 
defined macroscopically, is thus actually a property of a 
single spin, i.e. a strictly local property. 

As a second observable one can choose the occupa- 
tion probability, p, of the s z = level (averaged over all 
spins), 

P=-(El Q i><°il)- (37) 
i=i 

Similar to m z according to equation l|36(l , p may be writ- 
ten as 

P=(|Qfe)(0 fe |), (38) 
for any k — 1,2, ... ,n and is thus strictly local, too. 
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FIG. 10: m z (solid line) and p (dashed line) as a function of 
temperature T for a spin-1 chain of 4 particles. T is given in 
units of B and J = 0.1 x B. 




FIG. 11: m z (solid line) and p (dashed line) as a function of 
temperature T for a spin-1 chain of 4 particles. T is given in 
units of B and J = 2 x B. 



Now, if each single spin was in a canonical state with 
temperature 2| oc , m z and p would both have to be mono- 
tonic functions of T\ oc . In this case, T\ oc could, after 
calibration, be infered from measurements of m z or p, 
respectively. Note, that m z is proportional to the local 
energy, the average energy of one subsystem. 

Figure njj shows m z and p as a function of the global 
temperature T for a spin-1 chain of 4 particles with the 
Hamiltonian Q34|l for weak interactions, J = 0,1 X B, 
Both quantities are monotonic functions of each other. 

The situation changes drastically when the spins are 
strongly coupled. In this case the concept of temperature 
breaks down locally due to correlations of each single spin 
with its environment. 

Figure lTTl shows m z and p as a function of temperature 
T for a spin-1 chain of 4 particles with the Hamiltonian 
(|34JI for strong interactions J = 2 x B. Both quanti- 
ties are non-monotonic functions of T and therefore no 
mapping between m z and p exists. 

How could a local observer determine whether the sys- 
tem he observes, a single spin, is in a thermal state and 
can therefore be characterized by a temperature? The 
local observer would need to compare two situations: 

In the first situation, the spin is weakly coupled to a 
larger system, the heat bath. In this situation, the local 
observer could measure m z and p as functions of the tem- 



perature of the heat bath and would get a result similar 
to figure ^| This result would not be sensitive to the 
details of the coupling to the heat bath. The local ob- 
server would thus recognize this situation as a particular 
one and might term it the "thermal" situation. 

The second situation is fundamentally different. The 
spin is now strongly coupled to its surrounding. If the 
local observer again measures m z and p as functions of 
the temperature of the surrounding, he would get a result 
like in figure ITTI 

The observer can tell the difference between both situ- 
ations, even if he has no access to the global (true) tem- 
perature T of the surrounding. In the first case he can 
construct a mapping from say m z to p, i.e. p(m z ), or vice 
versa, m z (p), in the second he cannot: There exist for ex- 
ample two values of p corresponding to only one value of 
m z . Here the concept of a local temperature breaks down 
at least on the level of individual particles, since temper- 
ature measurements via different local observables would 
contradict each other. 

The question, whether and on what scale local temper- 
atures can exist in systems, which are in a global equi- 
librium state is thus indeed physically relevant. The ad- 
vantage of the concept of temperature is that it allows 
to predict various physical properties of the considered 
system. This is only possible if different properties (ex- 
pectation values of observables) map one to one on each 
other as in figure The following example illustrates 
the situation: 

Consider a piece of metal, say a wire. Assume , its tem- 
perature is measured via its electrical resistance. Why 
are we interested in this temperature? We are interested 
in it because it also allows us to predict how the wire 
behaves with respect to other physical processes. For ex- 
ample, if we know its temperature, we can tell whether 
the wire is going to melt or not. Effectively, we thus have 
a mapping of the resistance onto the fact that the wire 
is going to melt or is not going to melt. In a more math- 
ematical language, we can construct a function: melting 
as a function of the resistance. Analogously, for the sce- 
nario of figure a local observer is able to construct a 
function m z (p). 

What happens if such functions can no longer be con- 
structed? In this situation, the concept of temperature 
becomes useless. Assume our wire had such properties. 
We could still measure its resistance and, if we wished, 
could assign a "temperature value" to it. This "temper- 
ature value" , however, would be of no further use, since 
it would not allow us to predict whether the wire is go- 
ing to melt or not. A situation where such problems do 
really occur is the scenario of figure ITTI 



C. Potential Experimental Tests 

Finally, we address the question of whether the effects 
described here could be observed in real experiments. In- 
deed, pertinent experiments are available and have partly 
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already been carried out: 

A realization of a quasi one dimensional anti- 
ferromagnetic spin-1 Heisenberg chain is the com- 
pound CsNiCl 3 HI HI- Here the coupling is 
J 2.3 meV. To achieve a detectable modulation of 
m z and p, the spins should be significantly polarized for 
T > 0. Therefore a sufficiently strong applied magnetic 
field is needed. For CsNiCl3, a field of roughly 9.8 Tesla 
would correspond to J = 4 x B. 

The magnetization in an applied field can be measured 
with high precision by means of a SQUID j^. The occu- 
pation probability of the s z = states, on the other hand, 
is accessible via neutron scattering experiments j3^, HtJ . 
The differential cross section for neutron scattering of 
spin-systems is a function of the Fourier transforms of 
spin correlation functions |4(ij |. One can thus obtain in- 
formation about the quantity 

- J2wm°m+ ^(o)^(o)) = i+ P (39) 

r 

from the measurement data. Therefore, p is measurable 
in neutron scattering experiments. 

Such experiments or a combination thereof could thus 
be used to demonstrate the non-existence of local tem- 
perature. 

IX. CONCLUSION AND OUTLOOK 

In the present work, we have considered the mini- 
mal spatial length scales on which local temperature can 
meaningfully be defined. For large systems in a global 
equilibrium state, we have reviewed the derivation of 
two criteria which are valid for quantum many body sys- 
tems with nearest neighbor interactions and discussed the 
physical relevance of the existence and non-existence of 
local temperatures. 

Some questions related to the microscopic limit of the 
applicability of thermodynamics have thus been clarified. 
Nevertheless, some open problems remain and even new 
ones appeared in the context of the present approach. 

First of all, the generalization of the calculations to 
scenarios with only local but no global equilibrium is an 
issue of significant importance. One might expect that 
the length scales do no longer depend on the interactions 
and the global temperature only, but that the tempera- 
ture gradient becomes relevant, too. 

For global non-equilibrium, local temperature mea- 
surements of the standard type are very interesting and 
important issues on their own. As we have discussed 



here, these measurements have no spatial resolution if 
the sample is in a global equilibrium state. On the other 
hand, local temperature measurements with spatial reso- 
lution are being done for macroscopic setups and nobody 
would dare to question their validity. Therfore the max- 
imal spatial resolution of this kind of measurements is 
an interesting question and the present understanding of 
this topic is quite poor. 

Future research could also be concerned with new 
physics that might show up for small entities, which are 
in contact with a thermal surrounding, but show non- 
thermal behavior due to the breakdown of temperature 
on the respective scale. One example for this are the ob- 
servable features discussed in section IVIIII However, one 
might think about more surprising phenomena, as for ex- 
ample anomalous pressure fluctuations in very small gas 
bubbles enclosed in a piece of solid. With respect to fu- 
ture nanotechnologies, such phenomena could equally be 
harmful or useful, depending on whether one is able to 
design the devices in the pertinent way. 

Finally, possible generalizations of thermodynamics, 
that could apply on even smaller scales seem to be in- 
teresting. In the present work, we have considered the 
microscopic limit of usual thermal behavior in quantum 
systems, i.e. Quantum Thermodynamics |4l|. where ef- 
fective interactions among the considered parts are small. 
One might thus wonder, whether only partitions with 
weak effective couplings can be considered within such 
a "universal" description, that does not depend on the 
details of the microscopic constituents, or whether there 
exists again an intermediate level of description, not as 
universal as standard thermodynamics but applicable on 
smaller scales. Since, in standard thermodynamics, equi- 
librium states are fully characterized by one single pa- 
rameter, temperature (cf. equation ©), one could for 
example imagine that there exists a class of generalized 
equilibrium states, which require say two or three param- 
eters for their characterization. Some phenomenologi- 
cal attempts in this direction have already been made 
[42l l43j ]. Nonetheless, a justification of these attempts 
from an underlying theory, i.e. quantum or classical me- 
chanics, is still missing. 
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